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Abstract 

> ' 

. The poor man's scaling is extended to higher order by the use of 

the open-shell Rayleigh-Schrodinger perturbation theory. A generalized 
. Kondo-type model with the SU(n)xSU(m) symmetry is proposed and 

renormalized to the third order. It is shown that the model has both 
local Fermi-liquid and non-Fcrmi-liquid fixed points, and that the latter 
' becomes unstable in the special case of n = m = 2. Possible relevance of 

the model to the newly found phase IV in Ce s Lai_ x B 6 is discussed. 



1 Introduction 



In the renormalization theory one seeks a small parameter to control pertur- 
bative calculation of physical quantities. In the Kondo model, however, one 
encounters unmanageable divergence for physical quantities because the fixed 
point of the model corresponds to the strong coupling limit. Thus the numerical 
renormalization group has proved useful for quantitative study of the Kondo 
effect |TJ. On the other hand, the multi-channel Kondo model has a nontrivial 
fixed point at the exchange interaction J c which becomes small as the number 
of channels increases. Then the model has a controlled perturbation expansion 
from the limit of large degeneracy of conduction bands [§, ||) . It is highly desired 
to have a simple scheme to perform perturbative renormalization to any desired 
order. This paper proposes an extension of the poor man's scaling to higher or- 
ders by the use of the Rayleigh-Schrodinger perturbation theory together with 
the open-shell formalism. 

Another motivation of the present study is due to observation of a strange 
phase (called IV) in Ce^Lai-^Bg 0, [|, [£§. The magnetic susceptibility shows 
a cusp on entering the phase IV from the paramagnetic phase with decreasing 
temperature. This suggests that the Neel state is present here. In contrast to 
the phase III which has both antiferromagnetic and quadrupole orders, however, 
the phase IV has a very small magnetic anisotropy and almost no magneto 
resistance. We recognize the importance of the orbital Kondo effect which is 
active even in the presence of a spin ordering. In order to inspect the importance 
one needs to know how the intersite interaction and the on-site orbital and spin 
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Kondo effects compete. As the first step toward this direction, we study the 
simplest impurity model that displays both the orbital and spin Kondo effects 
with special attention to the combined spin and orbital symmetry of the 4/ 
shell. 

In conventional perturbation theory for Kondo-type models one introduces 
fictitious fermions or bosons to represent the localized spin |8], |9j. Then one 
computes some lower-order Feynman diagrams in order to apply the multiplica- 



tive renormalization [|1C], 11]. However, it is known that there is no linked-cluster 
expansion in this pseudo-particle theory [O]. Thus in higher-order it is diffi- 
cult to accomplish proper counting of perturbation processes. To overcome the 
difficulty several methods have been proposed. Among them is the resolvent 
formalism or the Brillouin-Wigner-type perturbation theory, which does not 
rely on the linked-cluster property [13, 14, 15, 16]. Another elegant way is to 
introduce a particular value of the imaginary chemical potential for the ficti- 
tious fermions. Then unphysical contributions cancel each other [17, 18]. We 
note that the latter method works only for the magnetic impurity with spin 
1/2- 

In doing perturbative renormalization of the Kondo model, the "poor man's 
scaling" of Anderson [19| provides the simplest framework for practical compu- 
tation. If one wants to proceed to higher-order renormalization, however, ex- 
tention of the scaling method is not straightforward |2(], 21]. In this paper, we 
propose an alternative approach to extend the poor man's scaling to arbitrary 
high orders. In contrast to the original idea of using the invariance property 
of the t-matrix, the present theory considers the effective Hamiltonian which 
gives the same energy spectrum as the original one at each order of renormal- 
ization. In deriving the effective Hamiltonian we use the Rayleigh-Schrodinger 
perturbation theory in the open-shell formalism, which is popular in quantum 
chemistry and nuclear physics p2| , |23|. We demonstrate the usefulness of our 
scheme by computing the third-order term in the renormalization group. 

The organization of the paper is as follows: In the next section we review 
the effective Hamiltonian formalism by the use of the Rayleigh-Schrodinger per- 
turbation theory. Section 3 performs renormalization to third order explicitly 
for the multi-channel Kondo model. We explain the notion of the folded di- 
agram which is very convenient to exploit the linked cluster property. These 
two sections do not give new results to the Kondo problem itself, but explains 
the formalism and demonstrates its simplicity. In Section 4 we introduce an 
SU(n)xSU(m) model where n denotes the number of spin degrees of freedom, 
and m does the orbital one. This model is different from the SU(n)xSU(m) 
Anderson model already discussed in the literature [^4|, pq] . Our model in- 
cludes both the Fermi-liquid and non-Fermi-liquid fixed points, and reduces to 
the multi-channel Kondo model in a special case. The renormalization to the 
third order is performed for general n and m by the use of the Lie algebra. We 
give detailed analysis of each fixed point and the renormalization flow. The 
final section summarizes the paper and describes possible further development. 



2 Effective Hamiltonian in the Model Space 
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2.1 Projection to the model space 

As a preliminary to perturbative renormalization we introduce the projection 
operator P to the model space M. We divide the Hamiltonian H = Hq + V 
such that the unperturbed part Hq commutes with P. In many cases we are 
interested only in certain lower part of the full spectrum of H. The effective 
Hamiltonian H e ff is designed to give the identical spectrum in the model space 
as that of H. Namely if ipi is an eigenfunction of H with the eigenvalue Ei, we 
require 

H eff P^i = EiP^i. (1) 

To construct H e ff it is necessary to introduce the wave operator f2 which re- 
produces tpi from Ptpi, namely tpi = fiP^. Note that the operator product QP 
is not an identity operator since P does not have an inverse. We emphasize 
that the above reproduction is possible only because ipi is an eigenfunction of 
H. 

In the Brillouin-Wigner perturbation theory the wave operator can be de- 
rived in the simple form p2|: 



n=0 



where Q = 1 — P is the projection operator to the space complementary to 
the model space. In spite of its simple appearance, the energy dependence of 
the wave operator makes it rather inconvenient for term-by-term perturbation 
theory in higher orders. On the other hand, the Rayleigh-Schrodinger pertur- 
bation theory is able to derive and hence H e ff i n a form which involves only 
eigenenergies of Hq and is free from the unknown energy E^. We write the 
Schrodinger equation in the two forms: 

(E - Ho)ij) = Eip- H ttPiJj = VnPip, (3) 

n(E - H )Pip = e^- nH Pip = npvnpip, (4) 

where the index i to specify an eigenstate is omitted. The latter equation is 
obtained by application of P on both sides of eq.(|3|) and further application of 
Q. Subtracting eq.(Q) from eq.(|3|), we obtain 

[Q, H }Pip = (1 - tlP)VnPifj. (5) 

We make the power series expansion: 

n = n + ni + n 2 + ■ ■ ■ , 

where Q n denotes the 0(V n ) contribution with Qq = 1. Comparing the terms 
with the same order of magnitudes on both sides of eq.(||), we obtain 

n-l 

[n n , H ] = QVn n -i - J2 tojPVSln-j-x. (6) 

3=1 

The effective Hamiltonian is accordingly expanded as 

H eff = P(H + V)P + H 2 + H 3 + ..., (7) 
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with H n = PVn n -iP for n > 2. 

Matrix elements of some lower-order terms are explicitly given by 

(a\H 2 \b) = (a\V{e b -H Q )- l QV\b), (8) 
(a\H 3 \b) = (a\V^—QV^—QV\b)-^2(a\V^—^—QV\c)(c\V 

Q> — -Ho £6 — -Ho c £fe — H e c — Hq 

where states such as \a), \b) and |c) belong to the model space with Ho\b) = eb\b} 
etc. We note that the third-order Hamiltonian is not Hermite. It is possible to 
convert it to a Hermitian operator by a similarity transformation |23[| . 

2.2 Linked cluster property 

In the Rayleigh-Schrodinger perturbation theory, two kinds of terms appear in 
H n with n > 3; the first kind takes the simple form [(e& — Hq)~ 1 QV] ti and the 
other kind of terms have less simple form and originate from the second term 
in eq.(||). The latter terms in fact play an important role in leading to the 
linked-cluster expansion. Namely, among all processes appearing as a result of 
Wick decomposition, unlinked diagrams are cancelled by the second kinds of 
terms |2^]. Here the linked or unlinked diagrams have the same meanings as 
those in the Feynman diagram method. 

The Rayleigh-Schrodinger theory has two cases to treat: the simpler one 
is the closed shell case where a large energy gap separates the occupied and 
unoccupied fermion orbitals. The other one is the open shell case where in 
addition to the stable core orbitals, a part of nearly degenerate orbitals are 
occupied. The latter orbitals are called valence orbitals, and constitute the 
model space. The essential point of applying the open-shell formalism to Kondo- 
type models is to regard the local electron as belonging to the set of valence 
orbitals. In order to utilize the linked-cluster property of the perturbation 
theory, we remove the restriction that the number of local electrons be unity 
at any stage of intermediate states. The final form of the effective Hamiltonian 
conserves the number of local electrons in each order of perturbation. Hence 
in the model space with one and only one local electron which obeys either 
the Fermi or Bose statistics, one can reproduce the actual situation in the 
Kondo-like impurity model with the local number constraint. Note that this 
reproduction is possible because we are working with the canonical ensemble 
instead of the grand canonical one. 

In order to deal with the energy denominators of terms such as the second 



term in eq.(||), the concept of the "folded diagram" is useful [22, 23 1. Namely 
each time the projection operator P appears in the expansion, one draws a 
diagram where the direction of propagation is reversed. Then the energy de- 
nominator is associated with difference of left-going energies and the right-going 
ones just as in the first term of eq.@. By comparing the energy denominators 
given by the second term of eq.(^), and those given by the rule for the folded 
diagrams, one can confirm that both give identical results. 
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3 Renormalization of the Mult i- Channel Kondo Model 



As the simplest application of the present method, we consider the multi- 
channel Kondo model defined by 

H K = T,t,"A la c kla + JS-s (10) 

with channel index I. Here S is the localized spin with magnitude 1/2 and 

1 n 
kk' <*P '=i 

with N being the number of fc-states and er the Pauli matrix. It is known that 
the fixed point of the model gives the non- Fermi liquid for n > 2 |§] . In order 
to perform renormalization explicitly we assume the constant density of states 
for n conduction bands each of which per spin is given by 

p c (e) = (2Dy l 

for |e| < D and zero otherwise. The model space consists of the local spin and 
such conduction electrons whose energy e is within the range [— D + \dD\,D — 
\dD\] where dD(< 0) is the infinitesimal change of the cut-off energy. In order to 
simplify the diagrams we define the vacuum as the Fermi sea of non-interacting 
conduction electrons, and work with the particle- hole picture for excitations. 

In the lowest order, we consider the two processes shown in Figjl]. We 
associate the energy ±D to the intermediate conduction states which belong to 
the orthogonal subspace projected by Q. As is clear from the direct calculation, 
the effective exchange interaction depends on the energy of conduction electrons 
in contrast to the bare interaction. We replace this energy dependent interaction 
by a constant, and continue the renormalization with the same form of the 
effective Hamiltonian. Namely it is assumed that the dimensionless effective 
interaction g{D) with a given cut-off energy D is related to another with a 
different cut-off energy D' by a multiplicative factor z(g,D' /D). We require 

g(D , ) = z(g,D'/D)g(D). (12) 

This multiplicative property can be expressed equivalently in the differential 
form: 

% = 0(g), (13) 

where I = InD and the beta- function in the right-hand side does not depend 
explicitly on I. This scheme is called the logarithmic approximation in the 



literature [26]. In the effective Hamiltonian approach, the logarithmic approxi- 
mation requires that the infinitesimal change of the effective interaction through 
one step of renormalization is proportional to D^ 1 and is independent of the 
external energy. In the lowest order the independence is approximately satisfied 
for the external energy much smaller than D. To proceed similarly in the third 
order, we set = in eq.(^) and deal with the effective Hamiltonian: 

Hi = V(H^QV) 2 + £ VH^(e c - H^QV^^V. (14) 
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(a) 



(b) 



Figure 1: Scattering processes in second order. The solid line shows a 
conduction-electron state, while the dashed line the local electron. The pro- 
jection operator Q requires the intermediate conduction-electron states to have 
energies near the band edges. 



Figure 2: Exchange scattering processes in the third order. The projection 
operator Q requires one of two conduction-electron states in the loop to have 
energies near the band edges. 

We refer to the literature p(| for detailed discussion of the logarithmic approx- 
imation by a different formalism. 

Before performing the third-order renormalization for the Kondo model, 
it is instructive to consider first the potential scattering from the impurity 
instead of the exchange scattering. If the particle-hole symmetry is present in 
the conduction band, we can generate a particle-hole conjugate diagram from 
any given diagram by reversing the direction of all the conduction-electron 
lines. Then the incoming conduction-electron line which corresponds to the 
annihilation operator in the effective Hamiltonian and the outgoing one for the 
creation operator are interchanged. The diagrams shown in Fig.[l] represent the 
simplest example of a pair of conjugate diagrams. For the potential scattering 
the two diagrams in Fig.[l] combine to zero. This is because an extra interchange 
of the incoming and outgoing fermion lines are involved in constructing the 
conjugate diagrams. To the contrary, the non-commutativity of spin operators 
lead to nonzero result from the sum of two diagrams in Fig|l|. 

In the third-order renormalization we first consider such contributions that 
become dominant for large n. These are diagrams which have a loop of conduc- 
tion electrons as shown in Figs.^ and |3], since each loop of conduction electron 
lines acquires the factor n by summing over degenerate orbitals. The diagram 
shown in Fig.|| corresponds to the first term in eq. (|14|) . The folded diagram 
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Figure 3: The third-order folded diagram. The assignment of energy denomi- 
nators are explained in the text. 

shown in Fig.|| corresponds to the wave-function renormalization in the Green 
function formalism. With the notion of folded diagrams one can confirm that 
the cancellation of potential scattering diagrams persists at least to the third 
order. Namely the diagram shown in Fig. 3] has the same magnitude but has 
the sign opposite to the one shown in Fig.|. Thus the two contributions cancel 
each other for the case of potential scattering. In both Figs.^ and|||, the energy 
denominators are obtained by associating the excitation energy D with one of 
two lines in the electron loop. The contribution is given by 

+ / 1 n ^ = 7v (15) 



D {-D + e 1 ) 2 Jo (-D-e) 2 D' 

The product of the spin operators takes the form S a S@ S" y Tv(s a s' y )s 13 where 
the trace is over the spin states of conduction-electrons. On the other hand 
the folded diagram has the same energy denominators as given by eq.([l5|), and 
the spin part is given by S a S' y S l3 Tv(s a s' y )s^ . Then contributions from the two 
diagrams combine to give 

Y / S a [S l3 ,S~ f }Tr(s a s r )s> 3 = ~S ■ s, (16) 

a/37 

where we use the identity [S a ,S^] = ie^^S" 7 with e a/ g 7 being the completely 
antisymmetric unit tensor. 

As will be shown shortly, all non-loop diagrams in the third order can be 
neglected for multiplicative renormalization. Thus we recover the known result 

I 

5 ■>'■■;■>■■ <"> 

where g = Jp c . A trivial fixed point of the model is g c = which, however, is 
unstable. Another fixed point of the renormalization group is given by g c = 2/n 
which is small in the case of large n. Thus the fixed point is within the reach of 
the perturbative renormalization. Linearization of the scaling equation around 
g c = 2/n shows that the latter fixed point is stable. 

Now we turn to non-loop diagrams in the third order. Figure ||] shows such 
contributions with incoming and outgoing conduction electron lines deleted. 
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Figure 4: Exchange scattering processes in the third order without a 
conduction-electron loop. The projection operator Q requires the double lines 
to have energies near the band edges. There are two ways of associating the 
direction of arrows for each diagram. 



These diagrams belong to the group called the "parquet" in the Green func- 
tion theory [26]. The two ways of attaching directions of arrows correspond to 
particle-hole conjugate diagrams. It can be seen by direct calculation that the 
diagrams shown in Fig.|](a) and (b) cancel each other not only for the poten- 
tial scattering but for the exchange scattering. This is because both have the 
same energy denominators and matrix elements but have different signs. On 
the other hand, the diagram shown in FigJ^(c) and (d) do not cancel among 
themselves. Moreover the diagram (d) is logarithmically divergent. This diver- 
gence is, however, not related to the Kondo effect since the potential scattering 
also gives this divergence. 

Let us ask the nature of these contributions to the effective interaction 
in the third order. They have a common feature of being constructed from 
the second-order diagrams through replacement of a bare interaction vertex by 
the second-order effective interaction. A diagram with such property is called 
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"reducible". We note that the change of the cut-off influences the shift of the 
ground-state energy, since the band-edge parts are no longer available to the 
effective Hamiltonian in the next step. To reproduce the identical energy shift as 
that of the original model, the effective Hamiltonian should contain a part which 
compensates the missing energy shift associated with the cut-off region. In the 
third order, this compensation comes from the expectation value of the reducible 
part of the effective interaction. In fact a part of the third-order energy shift 
comes from the diagram which is obtained from Fig.|| (d) by joining the outgoing 
conduction-electron line with the incoming one. The resultant energy shift 
corresponds to the expectation value of the logarithmically divergent effective 
interaction. In a similar manner expectation value of the diagram in Fig.|| (c) 
also compensates another missing part of the energy shift associated with the 
cut-off region. 

Thus we conclude that in the third-order non-loop diagrams do not con- 
tribute to the multiplicative renormalization. The effective interaction relevant 
to the renormalization group comes only from "irreducible diagrams" which 
cannot be constructed from lower-order diagrams through replacement of a 
bare interaction by a "vertex part" in the terminology of the Green function 
formalism. In the conventional renormalization theory, the number of diagrams 
increases considerably as one goes higher order in g. These terms have been 
discussed in the literature JO], |T^] with mutually inconsistent results. The in- 
consistency may partly be related to non-universal character involved in higher- 
order renormalization. We plan detailed account of 0(g ) renormalization and 
comparison with available exact results in a future publication. In this sec- 
tion we reproduced the old result just for the purpose of demonstrating the 
simplicity of the new formalism. 

4 Renormalization of Generalized Kondo-Type Mod- 
els 

4.1 SU(2)xSU(2) Kondo model 

The orbital dynamics in rare-earth systems such as Ce x L&i^ x T}Q has remarkable 
entanglement of magnetic and electric multipoles. The four-fold degenerate 
crystal field level T$ is often treated in the framework of the SU(4) symmetry. 
However, if one starts with the Anderson model and renormalizes off the charge 
fluctuation, the multiplet structure of 4/ 2 configurations act differently to the 
orbital and spin degrees of freedom || . Hence we consider the following model 

m- 

#2X2 = e k C kla C klo + H ex (18) 
kla 

where I = 1,2 denotes orbital index. The exchange part H ex is given by 

H ex = JpSf ■ Jsc + KpTf • eVc + AlJ2(PS^f)(ch^c). (19) 
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where S and s refer to spin matrices and T and r do the orbital ones. They 
are given explicitly by 

PSf = \Y.fl<r a pfl^ PTf=lT,fL*lrfra, (20) 

la/3 Ira 

f*S"Vf = \Y,Y,fl°>lrU (21) 

lr o/3 

for localized electrons. Similar notation has also been used for conduction 
electrons. There is a constraint to suppress the charge fluctuation of / electrons: 



£//a/*a = l- (22) 



la 

In the special case of J = K = /, the present model reduces to the Coqblin- 
Schrieffer model. To see this we note that the spin permutation operator V s 
and the orbital one Vi can be written as 

P s = 2/tS/-cWi (23) 

Pi = 2/+T/- C trc+i (24) 

Then the exchange part of the Coqblin-Schrieffer model is represented by 

I/t c t(2S. s + I)(2T.T + i)c/, (25) 

in an abbreviated notation. This interaction reduces to eq.(|T^) with J = K = I. 

Renormalization of the model proceeds in a similar manner as that in the 
Kondo model. In order to simplify the notation we use the convention that J, K 
and / denote the effective interaction and the unit of energy is such that p c = 1. 
In the second order cross terms of J and K vanish because the interactions with 
them commute each other. This is not the case with the I term. In the third 
order the relevant diagrams are again only those shown in Figs.|| and ||. The 
vertex parts there can be either J, K or /. On the other hand the energy 
denominators are the same as in the case of the Kondo model. The resultant 
scaling equations are given by 

°_± = _ (1 _J )( j2 + 3/ 2 )i (26) 
r)K 

— = -(1-K)(K 2 + 3I 2 ), (27) 
dl 

— = -2I{J + K)+I(K 2 + J 2 + 2I 2 ). (28) 
ol 

Let us discuss implication of the scaling equations. The set of equations has 
six fixed points (i) -(iv) as follows: 

(i) J c = K c = I c = 0, (ii) J C = I C = 0,K C = 1, (hi) J c = l,K e = I e = 0, 
(iv) J C = K C = 1, I c = 0, (v) J C = K C = 1, I c = 1, (vi) J C = K C = 1, I c = -1. 

The fixed points (ii), (hi) and (iv) correspond to the non-trivial fixed point 
known for the multi-channel Kondo model S which are in fact unstable in the 
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presence of I. We discuss in more detail the stability of these fixed points in the 
next section taking a generalized model. It is known that the Coqblin-Schrieffer 
limit of the model does not have the non-Fermi liquid. Then the fixed point (v) 
with I c = 1 should be an artifact of the third-order scaling. The correct fixed 
point is at I c = J c = K c = oo and gives the local Fermi liquid. 

A remarkable property of the scaling given by eq.(p8|) is the absence of 
the I 2 term. The SU(2) symmetry plays a special role in this absence as will 
become clear in the next section. Since the absence means that the sign of / 
is irrelevant in the present system, it suggests that the fixed point (vi) with 
I = —1 also flows to infinity in the exact renormalization group. This flow has 
been confirmed by the numerical renormalization group method |28|| . Namely 
the spectrum of the model with positive / is the same as that with negative I 
for the same absolute value. The flow to the Fermi liquid is also concluded by 



Pang [27], who started the renormalization from 1 = 0. 



4.2 SU(n)xSU(m) model 

Motivated by the previous success of the large n theory for the multi-channel 
Kondo model [||, ||, we now generalize the model of eq. (|l8|) to the SU(n)xSU(m) 
symmetry with arbitrary n and m. Before presenting the SU(n) xSU(m) model 
and its renormalization, we quote necessary formula for the Lie algebra of the 
unitary group SU(n). Let X a (A = 1, 2, . . . n 2 — 1) represent the set of generators 
of the SU(n) Lie algebra. The commutation rule is given by 

[X a ,X?] = if a ^X\ (29) 

where f a ^ is called the structure constant, and is completely antisymmetric 
against interchange of a pair of indices. In the SU(2) case f a p y reduces to the 
unit tensor e a/ g 7 of the third rank. Then we require the orthonormality 

Tr(X a X?) = ~5 a(3 , (30) 

which is the generalization of the relation obeyed by the SU(2) spin matrices. 
The completeness relation is expressed as 

Y;{X a U{X a ) cd = \(5adhc ~ -tabbed)- (31) 

a 2 n 

In the special case of SU(2), the above relation is translated to that of spin 
permutation V s of local and conduction electrons as given by eq. (B3h . 



By analogy with eq.(lS) we now introduce the model with n-fold degenerate 
spin degrees of freedom, and m-fold degenerate orbital degrees of freedom. The 
exchange part of the Hamiltonian is given by 

H ex = J^(/t^/)(ct^c) + ^^(/tT^/)(cVc) + 4/^(/t^T l '/)( C t s ^c). 

(32) 

where 

fa"/ = ~Y,fU X "UfiP, f^f=\Y,fl{X V )lrUa, (33) 

la(3 Ira 

ftSPVf = lY,Y,fL( Xfl U(X V )lrfrp (34) 

Ir a(3 
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for localized electrons. There is the constraint given by eq.fl22|) as in the case 
of the SU(2) x SU(2) system. 

In the special case of / = mJ/2 = nK/2, the present model reduces to the 
SU(n x m) Coqblin-Schrieffer model. This becomes apparent when the product 
VsVi is written in terms of generators of the Lie algebra. Namely we can write 
V s by using eq.pll) as 



P s = 2^(/tS"/)( c t/c) + -, 
, n 



(35) 



and similar generalization of eq. (|24|) . 

The scaling equation up to the third order is given by 



dJ 
dl 
dK 
~dl 

di 

01 



n 



1 



in 



-J 



J + 4 1 



1 



V 



2 V 2 



K 2 + 4 f 1 - 4 ) I 2 

n z 



-I(nJ + mK) 



mn 



2m 2n\ 2 

mn / 

n m J 

j 2 + k 2 + a[i-^-^]i 2 



(36) 
(37) 

(38) 



We note that eq.(]3^) has a finite I 2 term in general. The coefficient vanishes 
only if n = m = 2 with n and m integers. 

There are six fixed points of the scaling characterized by eqs.(|3"6|) - (|38[) . In 
correspondence to the SU(2) x SU(2) model, they are given by 



(i) Jc = 
(iii) J c 

(v) J r -- 



K C = I C = 0, 



= 2/m, K c 
2/m, K c 



0, 



2/n,I c 



(ii) Jc 
(iv) J c 
(vi) J c 



I c = 0,K c 
~- 2/m, K c = 
-- 2/m, K c = 



= 2/n, 
2/n,I c 
2/n,I c 



with 



m + n 



m 2 n 2 



m* — n* 

the absolute magnitude of which is much smaller than unity for m, n ^> 1. Ex- 
cept for the fixed point (vi), the perturbative renormalization is well controlled 
because all J C ,K C and |/ c | are small as compared with unity. The case (vi) be- 
longs to the strong-coupling regime and the third-order scaling is unreliable. In 
the Coqblin-Schrieffer limit, this fixed point actually flows to the Fermi liquid 
one in the exact renormalization group. 

Let us study the stability of the fixed points. The trivial fixed point (i) 
is unstable as in the case of the standard Kondo model. The fixed point (ii) 
constitute a saddle point in the K — I plane; it attracts the renormalization 
flow parallel to the K-axis, but repels it parallel to the /-axis. Similar is the 
case (iii) with J and K interchanged. The linearized scaling equations around 
the fixed points (iv), (v) and (vi) all contain a part: 



dJ_ 

dl 
dK 



nm ( J 

m 



nm I K 



1 (, 1 

— + (1 2 

m L 
1 



71- 



+ 1 



m' 



n 



, r 

2 I c 



(39) 
(40) 
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From this set of equations it turns out that the coefficient of the linear term 
in the right hand side is positive. This means that the fixed point attracts the 
renormalization flow parallel to either J- or i^-axes, independent of the value 
of I c . The flow parallel to the I-axis depends on I c as explained below. 

At the fixed point (iv) with I c = 0, the remaining linearized scaling equation 
reads 

31 (m + -) I. (41) 



dl \n m, 

The fixed point is unstable since the coefficient on the right-hand side is positive. 

On the other hand, the linearized scaling equation at the fixed point (v) 
with I c = I mn is given by 

dl mn(m 2 + n 2 ) 



dl m 2 n 2 — m 2 — n 2 



(I -Imn)- (42) 



Since the coefficient on the right-hand side is positive, the fixed point is stable. 
The fixed point describes a non-Fermi-liquid state. Within the third order scal- 
ing, this fixed point connects smoothly, as the degeneracies n and m decrease, 



to the one at J = K = —I = 1 which was derived in §4.1. However, the special 



symmetry in the case of SU(2) x SU(2) rejects this smooth connection in the 



exact theory [28]. Thus we suggest that the non- Fermi- liquid fixed point is 
stable except for the case of SU(2)x SU(2). It should be an interesting future 
problem to see how this strange behavior is related to the symmetry. 

Figure [| shows qualitatively the flow diagram of the renormalization group 
in the special case of n = m (3> 1) where the condition J = K , if present 
initially, continues to hold in all steps of renormalization. The finite effective 
interaction in the strong-coupling fixed point at I c = 1 has been replaced by 
the correct flow toward the local Fermi-liquid. 



5 Summary and Outlook 

In this paper we have proposed a new method to perform renormalization for 
Kondo-type models. In contrast to the standard approach which considers 
the vertex part, our method is concerned with the effective Hamiltonian. Our 
approach does not involve infinite-order summation, but still employ the log- 
arithmic approximation. In this way the multiplicative property in the renor- 
malization is assumed. 

In the SU(n)xSU(m) model nearly all the fixed points, except a strong- 
coupling one with J c = 1, are shown to be within the reach of the perturbative 
renormalization. One of the two non-Fermi-liquid fixed points, which has I c = 
I nm and is stable for large n and m, becomes unstable against the Fermi liquid 
one in the special case of the SU(2)xSU(2) symmetry. This is due to the 
unexpected symmetry between / and —I. Meaning of this hidden symmetry 
should be identified in a future study. 

The other problem to be pursued further is the renormalization in still 
higher order. In particular comparison with exact results obtained by the Bethe 
Ansatz [29, EH] and the conformal field theory [31] should be made. The latter 



approach to the multi-channel Kondo model derived the critical value of the 
exchange interaction as J c = 2/{m + 2), which tends for large m to the value 
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Figure 5: The renormalization flow of the SU(n)xSU(m) model in the special 
case of n = m. The asterisk represents the (fictitious) fixed point given by the 
third-order scaling. 

2/m obtained by the third-order renormalization. Since the critical value itself 
is not a universal quantity, the 1/m expansion does not necessarily reproduce 
the value J c = 2/{m + 2). However, observable quantities such as the critical 
exponent must be independent of the method of derivation. The Rayleigh- 
Schrodinger theory seems to be the simplest approach to compute the higher- 
order beta function. The details of the higher-order renormalization will be 
investigated in a separate paper. 

Concerning the physical quantities in the intermediate regime toward the 
fixed point, the spin susceptibility can have different temperature dependence 
from that of the orbital susceptibility. Since both spin and orbital fluctuations 
contribute to the resistivity, conventional argument to relate the transport and 
magnetic properties via the single Kondo temperature has to be revised in the 
presence of an orbital degeneracy. Recognition of this difference seems to a 
clue to understanding the nature of the phase IV of Ce^Lai-^Be. As another 
example of entangled spin and orbital degrees of freedom we mention transition 
metal oxides, especially manganites. For these systems we have to modify the 
model to account for the situation where the crystal field splitting is much larger 
than the spin-orbit splitting. We are now analyzing the SU(2)xSU(2) model 
quantitatively with use of the numerical renormalization group |2^] . 
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